® 



J) 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



©4Slumero de publication: 0 624 847 A1 



® 



© Numero de depot: 94201267.5 
@ Date de depot: 06.05.94 



DEMANDE DE BREVET EUROPEEN 

© int.ci.5 G06F 15/353, G06F 15/36 



@ Priorite: 12.05.93 FR 9305715 

0 Date de publication de la demande: 
17.11.94 Bulletin 94/46 

0 Etats contractants designes; 
DE FR GB IT 

© Demandeur: LABORATOIRES 
D'ELECTRONIQUE PHILIPS 
22, Avenue Descartes 
F-94450 Limeil-Brevannes (FR) 

@ FR 

© Demandeur: PHILIPS ELECTRONICS N.V. 



Groenewoudsew^eg 1 
NL-5621 BA Eindhoven (NL) 
© DE GB IT 

@ Inventeur: Deville, Yannick 
Societe Civile S.P.I.D. 
156, Boulevard Haussmann, F-75008 Paris 
(FR) 

© Mandataire: Landousy, Christian et al 
Societe Civile S.P.I.D. 
156, Boulevard Haussmann 
F-75008 Paris (FR) 



© Dispositif et methode pour generer une fonction a approximation. 
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© Dispositif (5) pour generer une fonction d'approximation fondee sur des prenniers couples ((X-.. Yi) a (Xe. 
Ys)) de valeurs associant une grandeur dependante (Y, a Ye) a une grandeur independante (Xi a Xe). et pour 
determiner des seconds couples (Xa, Y'a) de valeurs desdites grandeurs d'apres ladite fonction d'approximation. 
Le dispositif comprend : 

- des premiers moyens (10) pour determiner iterativemant au moins une fonction lineaire courante de 
regression, pour selectionner celle des fonctions lineaires courantes qui delivre ['approximation de tous les 
couples de ladite suite avec des erreurs minimales, et pour coder, a I'aide de codes specifiques (p. q). la 
fonction lineaire de regression selectionnee. 

- et des seconds moyens (17) pour determiner lesdits seconds couples (X^. Y'a) a I'aide desdits codes 
specifiques. 

Le dispositif peut etre utilise pour le calcul de valeurs approchees de fonctions mathematiques, par exemple 
dans un reseau de neurones, ou pour determiner une fonction de regression formant une approximation de 
valeurs experimentales de mesure, par exemple des mesures dispersees issues du suivi d'un processus 
industriel. 

L'invention concerne egatement une methode pour generer une fonction d'approximation. 
Application : calcul approche de fonctions mathematiques. calcul de fonctions de regression. 
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L'invention concerne un dispositif et une methode pour generer une fonction d'approximation. cetle-ci 
etant fondee sur des premiers couples de valeurs associant une grandeur dependante a une grandeur 
independante. et pour determiner des seconds couples de valeurs desdites grandeurs a partu de ladite 
fonction d'approximation. 

5 Un dispositif et une methode de ce type sont connus du brevet US-A- 3 789 203 qui decrit un 

generateur de fonctions operant une approximation par interpolation iterative. Ce dispositif est prevu pour 
etre utilise dans des applications de traitement de donnees necessitant un calcul de fonctions telles que 
sin(x), tg(x) par exemple. Ce dispositif ne requiert qu'une capacite de stockage minimaie de la part d'un 
dispositif utilisateur. A partir de deux points appartenant a une fonction a interpoler. la methode tout d'abord 
;o interpole la fonction par une droite reliant las deux points, puis fait une approximation des ecarts entre ia 
droite et la fonction par des approximations polynomiales de degre croissant. Ensuite elle substitue aux 
points initiaux des points approximatifs pour reduire ta longueur du segment reliant les points a traiter et 
enfin reitere les operations precedentes. 

Une telle methode necessite des ressources importantes en moyens de calcul et ne peut etre mise en 
75 oeuvre qu'avec des calculateurs performants. 

Or, il existe des applications ou une tel'e methode ne peut pas etre mise en oeuvre car elies ne 
disposent pas des ressources suffisantes. De plus, pour certaines applications on peut se satisfaire d'un 
calcul approche de ia fonction pour des valeurs, en nombre limite, de ta grandeur independante. 

II peut s'agir d'une fonction sigmoVde appliquee a des potentiels neuronaux delivres par au moins un 
20 neurone dans un reseau de neurones. II peut s'agir d'une autre fonction non lineaire. par exemple une 
fonction racine. pour des calculs de distances entre des etats de neurones. Les applications peuvent aussi 
concerner d'autres dispositifs comme des generateurs de fonctions, des calculateurs ou autres. 

Pour calculer une telle fonction. sans passer par une fonction d'approximation, on peut utiliser 
differentes manieres. 

25 On peut effectuer le calcul mathematique exact pour chaque valeur de la grandeur independante a 
traiter, en programmant un calculateur selon les methodes connues. Une telle methode necessite d'effec- 
tuer a chaque fois les memes operations ce qui peut necessiter beaucoup de temps si le nombre de 
valeurs est eleve. 

On peut aussi prealablement stocker dans une memoire des tables precalculees. Dans ce cas. la 
30 lecture en memoire du resultat peut etre rapide. fVlais pour couvrir avec un pas assez fin toutes les valeurs 
possibles de la grandeur independante. il faut alors disposer de tables de grandes capacites. Ces methodes 
de calcul presentent done des inconvenients. 

D'autre part, on peut etre conduit a identifier deux grandeurs qui sont dans la dependance I'une de 
Tautre par des couples de valeurs associant une grandeur dependante a une grandeur independante. Ainsi, 
35 dans le suivi d'un processus industriel, on peut etre conduit a mesurer par exemple un rendement R d'une 
operation en fonction de la temperature T a laquelle a ete realisee ladite opeation R = f(T). Un suivi du 
processus peut consister a consigner des lots de couples de mesures sur un graphique. Celui-ci peut etre 
utilise pour caracteriser le processus ou pour en deduire de nouveaux parametres de reglage de ladite 
operation. Ceci fait par exemple I'objet de I'article de H. ISHIBUCHI et H. TANAKA, "Regression analysis 
40 with interval model by neural networks" in "IEEE International Joint Conference on Neural Networks", vol. 2. 
18-21 Nov. 1991, SINGAPORE. II est souhaitable que ces nouveaux parametres soient representatifs du 
fondement de ladite operation et que les fluctuations de mesures inherentes a ce genre de processus 
soient ecartees. II est done souhaitable de determiner une approximation de la fonction f(.). 

Ainsi dans un cas il peut s'agir de mesures erratiques ou entachees d'erreurs que Ton desire 
45 representor par une fonction d'approximation. 

Dans un autre cas, on connait des valeurs precises mais I'utilisation a en faire ne necessite pas une 
grande precision et une fonction d'approximation suffit. 

Un des buts de invention est de generer une fonction d'approximation avec des moyens materiels 
reduits permettant de calculer rapidement un nombre reduit de valeurs de la grandeur dependante utiles a 
50 I'application sans avoir pour cela a determiner d'autres valeurs de la fonction d'approximation. Un but 
complementaire est de delivrer des valeurs qui peuvent etre approchees dans la limite d'une erreur 
maximale controlee. 

Ce but est atteint avec un dispositif caracterise en ce qu'il comprend : 
- des premiers moyens : 

55 - pour determiner iterativemant au moins une fonction lineaire courante de regression en rendant 

egales, en valeur absolue. des premieres erreurs de signes alternes mesurees entre des premieres 
valeurs de la grandeur dependante pour trots couples d'une suite desdits premiers couples, et 
respectivement des secondes valeurs de la grandeur dependante determinees, d'apres ladite 
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fonction lineaire courante. pour les memes valeurs de la grandeur independante. 

- pour selectionner celle des fonctions tineaires courantes qui delivre rapproxinnatton de tous les 
couples de ladite suite avec des erreurs minimales. 

- et pour coder, a I'aide de codes specifiques. la fonction lineaire de regression selectionnee, 

5 - et des seconds nnoyens pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 

AinsI avantageusement on determine une fonction lineaire de regression approchant au mieux les 
differents couples de valeurs connus. Les resultats approches ainsi delivres forment un compromis 
satisfaisant pour de nombreuses utilisations du dispositif generateur de fonctio... 

Une fonction lineaire de regression est une fonction simplificatrice qui represente un phenomene 
w complexe en reduisant les parametres significatifs. En representant la suite de couples de valeurs par des 
points dans un espace a deux dimensions, la fonction lineaire de regression devient une droite de 
regression. 

Ainsi apres avoir defini la droite de regression par des codes, on peut calculer une valeur approchee de 
la grandeur dependante en tout point de la droite de regression avec des moyens reduits pour des valeurs 
75 quelconques de la grandeur independante. 

L'invention concerne egalement une methode pour generer une fonction d'approximation, la meihode 
comprenant : 

- une premiere phase : 

- pour determiner iterativement au moins une fonction lineaire courante de regression en rendant 
20 egates, en valeur absolue. des premieres erreurs. de signes alternes. mesurees entre des 

premieres valeurs de la grandeur dependante pour trois couples d'une suite desdits premiers 
couples, et respectivement des secondes valeurs de la grandeur dependante determinees, d'apres 
ladite fonction lineaire. pour les memes valeurs de la grandeur independante. 
pour selectionner cede des fonctions lineaires courantes qui delivre I'approximation de tous les 
25 couples de ladite suite avec des erreurs minimales, 

- et pour coder la fonction lineaire de regression selectionnee a I'aide de codes specifiques. 

- et une seconde phase pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 
Les moyens mis en oeuvre par l'invention peuvent etre formes par un calculateur programme ou par un 

circuit dedie. lis peuvent aussi mettre en oeuvre des neurones. 

30 Un dispositif mettant en oeuvre des neurones selon l'invention peut etre utilise par un reseau de 
neurones, dont il peut notamment en constituer un sous-ensemble. En effet, pour foncttonner. le reseau de 
neurones doit disposer de moyens pour appliquer une fonction non lineaire d'activation aux potentiels de 
neurones qu'il delivre. Selon l'invention, le dispositif muni de neurones peut calculer une approximation de 
cette fonction non-lineaire d'activation. II peut egalement calculer des distances entre des etats de neurones 

35 en catculant une approximation d'une fonction racine carree destinee a etre exploitee dans le reseau de 
neurones. 

Lorsque la taille de la suite de couples de valeurs fournis initialement est elevee, on peut diviser la suite 
de couples en plusieurs sous-ensembles pour determiner plusieurs droites de regression et ameliorer la 
precision de I'approximation. La fonction d'approximation de la suite de couples est alors formee par une 
40 fonction lineaire par morceaux pour laquelle une exigeance de continuite entre les morceaux peut etre ou 
non imposee. 

Certains couples de I'ensemble de couples de valeurs peuvent avoir une influence particuliere que Ton 
peut concretiser en donnant un coefficient de ponderation specifique a chaque couple. Dans ce cas. I'erreur 
affectee a chaque couple tient compte de ce coefficient de ponderation specifique. 
45 Ces differents aspects de l'invention et d'autres encore seront apparents et elucides a partir des modes 
de realisation decrits ci-apres. 

L'invention sera mieux comprise a I'aide des figures suivantes donnees a titre d'exemples non limitatifs 
qui representent : 

Figure 1 : un graphique montrant une representation a deux dimensions d'un ensemble de points avec 
50 une droite de regression D. 

Figure 2 : un graphique montrant un ensemble de points et des droites servant a la determination d'une 
enveloppe. 

Figure 3 : un organigramme d'une premiere variante de mise en oeuvre de la methode a partir de 
triplets de points. 

55 Figure 4 : un organigramme d'une seconde variante de mise en oeuvre de la methode a partir de triplets 
de points. 

Figure 5 : une partie d'organigramme d'une troisieme variante de mise en oeuvre de la methode a partir 
de couples de points. 
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Figure 6 : une partie d* organigramme indtquant une preselection de points appartenant a une enveioppe 
inferieure ou a une enveioppe superieure de I'ensennble des points. 

Figure 7 : un organigrannme pour la determination des enveloppes infeneure et superieure. 
Figure 8 : un schenna d'un dispositif selon Tinvention. 
5 Figure 9 : un schema d'un dispositif de codage pour le calcul et le codage d'une droite de regression 
dans le cas general. 

Figure 10 : un schema d'un dispositif de codage pour le calcul et le codage d'une droite de regression 
pour la premiere variante. 

Figure 11 : un schema d'un dispositif de transcodage pour le calcul des valeurs de la grandeur 
10 dependante a partir d'un ensemble de droites de regression codees. 

Figure 12 : un schema d'une realisation neuronale des moyens de comparaison de la figure 10. 
Figure 13 : une representation montrant une fonction d'approximation formee de plusieurs droites. 
Figure 14 : deux representations montrant une determination des limites de validite de deux droites 
consecutives. 

15 Figure 15 : une representation concernant le raccordement de droites consecutives. 

L'invention concerne I'approximation d'une fonction connue uniquement a travers un certain nombre de 
points par exemple Pi... Ps (figure 1) dans une representation a deux dimensions. Chacun de ces points 
est defini par un couple de valeurs (x. y) reliant la grandeur independante x a la grandeur dependante y. 
Par la suite, il sera question d'une paire (respectivement d'un triplet) de points, ce qui met en oeuvre deux 
20 couples (respectivement trois couples) de valeurs. On ordonne les points d'apres un ordre croissant des 
valeurs d'abscisses Xi, ce qui definit un indice i croissant avec lesdites valeurs. Une convention inverse 
peut etre faite en modifiant en consequence I'expose ci-apres. 

Selon Tinvention on effectue une approximation de I'ensemble des couples (Xi, Yi), (X2. Y?),.. par une 
droite de regression D ayant pour equation : 

.25 

(1) D : y = p.x + q, 

ou x et y sont des variables courantes. 

Pour cela on considere dans une premiere variante trois couples de valeurs, par exemple (X3. Y3), (Xa. 

30 Y4). (Xs. Ys). et on determine une droite de regression D par equilibrage des ^^reurs absolues. Une erreur 
est mesuree par la difference apparaissant, pour une abscisse x donnee. entre la valeur y du point ei 
Tabscisse y mesuree sur la droite de regression. Equilibrer tes erreurs sur trois points, consiste a avoir trois 
erreurs egales en valeur absolue avec un signe d'erreur oppose aux deux autres pour le point ayant une 
abscisse x comprise entre les abscisses x des deux autres points. Puis on examine si pour les points 

35 restants de I'ensemble. I'erreur qui tes separe de la droite reste inferieure ou egale en valeur absolue a 
I'erreur prealablement determinee pour les trois points selectionnes. Ceci decoule du fait que Ton 
s'interesse a une erreur pire cas relative a I'ensemble de tous les points a prendre en consideration c*est- 
a-dire la plus grande erreur. en valeur absolue. qui existe entre un des points et la droite de regression. Si 
toutes les erreurs sont effectivement inferteures ou egales. la droite est selectionnee pour representer les 

40 points sinon on recommence les operations avec trois nouveaux couples de valeurs pour determiner une 
autre droite de regression. 

II peut exister plusieurs droites de regression representant tous les points de I'ensemble. Selon la 
methode on determine la droite de regression optimale qui minimise I'erreur pire cas definie prealablement. 
La figure 1 represente un exemple forme de six points Pi a Pg disposes selon une representation a 

45 deux dimensions. A des fins d'explication considerons le resultat final. On observe que la droite de 
regression D de la figure 1 est situee de sorte que les erreurs sont egales en valeur absolue pour les points 
P3, P* et Ps. Pour les points Pi. P2. Ps les erreurs sont inferieures aux precedentes. en valeur absolue. 
Dans un ensemble de points Pi a Pb. la methode va ainsi consister a rechercher les trois points 
particuliers. ici P3. P4. Ps. qui permettent de determiner la droite de regression optimale qui minimise 

50 I'erreur pire cas puis a coder cette droite. Dans le cas de la figure 1 representant un resultat final, si on 
trace deux droites Di et D2, paralleles a la droite de regression D, qui passent respectivement par les 
points P3. Ps. d'une part et Pd d'autre part, on constate que tous les points de I'ensemble sont a I'interieur 
d'un bandeau limite par les droites Di et D2 ou sur ces droites. 

La phase de determination de la droite de regression peut donner lieu a plusieurs mises en oeuvre dont 

55 seules les plus avantageuses seront decrites ci-apres. 
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Premiere variante de la premiere phase de la methode . 

La figure 3 represente la suite des etapes a mettre en oeuvre pour determiner la droite de regression 
optimale. 

5 Parmi Tensemble des points Pi... Pu (bloc 100). on selectionne (bloc 102) trois points quelconques P., 
Pj. Pk avec i < j < k. Ces trois points servent a determiner la droite de regression qui minimise I'erreur sur 
y pour ces trois points. Cette determination est effectuee de maniere analytique preferentiellement par des 
moyens programmes. On determine une droite de reg.ession D pour que trois erreurs (bloc 104) : 
EpD (P;. D). EpD (Pj, D). Epo (Pk. D) 

10 entre la droite D et chacun des points verifient : 

Epo (P.. D) = -EpD (Pj. D) = EpD (P..D) 



;5 



20 



25 



30 



avec 



Epo (Pi. D) = Yi - (p.Xi + q) 

et des relations analogues pour les autres erreurs. 

On determine la droite de regression D a I'aide des coefficients p et q de Tequation (1) tels que 



P = 



^k - 
^k - 



et 



2(X^ - X,) 



L'erreur se rapportant a un triplet (Pi, Pj. P^) s'ecrit alors : 
35 Er (P, P,, P,) = I Epo {P, D) I 

Lorsque l'erreur Et a ete ainsi calculee, on examine si les autres points de ('ensemble engendrent des 
erreurs inferieures ou egales, en vaJeur absolue. a celles des points Pj. Pj, P^. Pour cela on selectionne un 
point Prry (bloc 106) additionnei et on calcule la valeur absolue de l'erreur Epm (bloc 108) apparaissant entre 

40 la valeur de la variable y au point P^ et la droite D. 

Lorsque cette erreur Ep,^* est inferieure ou egale en valeur absolue a Et (bloc 110) (repere Y). le point 
Pm additionnei est accepte et la methode se poursuit (bloc 112) avec un point additionnei suivant (bloc 
106). Si tous les points additionnels satisfont au critere | Ep^ | ^ Er . ta droite D est acceptee et ses 
coefficients sont utilises pour coder la droite de regression optimale Dop, (bloc 114). 

45 Lorsque cette erreur Ep„, est superieure a Et (bloc 110) (repere N), le triplet de points Pj. Pj, P^ 
selectionne n'est pas accepte et un autre triplet de points (bloc 116) est choisi (lien 101) dans I'ensemble 
de points (bloc 102). La methode s'acheve a t'obtention de la droite D (bloc 114) satisfaisant ce critere 
meme si tous les triplets de points n'ont pas ete examines, 

II peut apparaitre des situations ou a Tissue de I'etape 116 tous les points possibles ont ete examines 

50 et aucun des triplets n'a fourni de solution (bloc 118). Dans ce cas, il est possible de reprendre le 
deroulement de la premiere variante en remplagant le test du bloc 1 10 par le test suivant : 

I Ep^\<a£r 

55 ou Q est un coefficient legerement superieur a 1. Dans cette meme situation, il est egalement possible de 
faire appet a la deuxieme variante. 

La premiere phase peut comprendre les etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite. 
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B - calcul de la foncnon Imeaire courante de regression D et determination d'une erreur de triplet 

Et-\Epd^ s'y rapportant. 
C - selection d'un couple additionnel. 

D - catcul d'une erreur additionnelle Epm entre le couple additionnel et ladite fonction. 
5 E - et lorsque | Epm \ ^ Et{\'\0) pour le couple additionnel. la methode reprend a I'etape C avec un 
couple additionnel suivant, 
F - et lorsque | £pm | > Et pour au moins un couple additionnel. 
la methode reprend a i'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite 
de couples, 

w G - et lorsque | Ep^ \ ^ Et pour tous les couples additionnels, la fonction lineaire courante de 
regression est codee et stockee en tant que fonction lineaire d'approximation. 
On peut choisir de scruter I'ensemble des triplets en prenant un ordre croissant ou un ordre decroissant 
ou un ordre ateatoire pour effectuer cet examen. Le triplet qui sera retenu pour determiner la droite de 
regression pourra de ce fait etre detecte a un instant quelconque du deroulement de cette scrutation. 11 

75 s'ensuit que la rapidite d'obtention de la droite de regression depend de I'instant au cours duquel le triplet 
est detecte. Sa mise en oeuvre presente un degre de complexite allant de N a N* ou N est le nombre de 
points initiaux. Sa complexite est done reduite pour un petit nombre de points. Cette variante permet 
d'obtenir une realisation matertelle avec une forte paralleiisation. Elle est tres peu sensible a une troncature 
des valours et fournit un resultat exact. 

20 

Deuxieme variante de la premiere phase de la methode. 

Dans cette deuxieme variante (figure 4). on selectionne successivement des triplets de points, on 
calcule chaque fois une droite de regression, et, par recurrence, on selectionne celle qui delivre I'erreur Epo 
25 la plus grande c'est-a-dire correspondant a I'erreur pire cas pour I'ensemble de points consideres. 
Selon la seconde variante, la premiere phase comprend tes etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite. 

B - calcul d'une fonction lineaire courante de regression D et determination d'une erreur de triplet 
Et=\Epo\ s'y rapportant, 

30 C - comparaison de I'erreur Et avec une erreur optimate Eop initiaiisee a une valeur strictement 
negative. 

D - et lorsque Et > Eop . mise a jour de I'erreur optimate Eop en remplagant Eop par Et et mise a jour 
des codes d'une fonction lineaire optimale de regression Dop en remplagant ceux-ci par les 
codes de la fonction lineaire courante D, 
35 E - puis retour a I'etape A pour selectionner trois autres couples, 

F - et lorsque tous les triplets de couples de valeurs de la suite ont ete testes, les derniers codes de 

la fonction lineaire optimale Dop constituent les codes de la fonction lineaire d'approximation. 
A chaque examen d'un triplet, on compare I'erreur Et dudit triplet avec I'erreur optimale precedemment 
memorisee et on met a jour I'erreur optimale Eop avec ia plus grande valeur de I'erreur Et determtnee pour 
40 chaque triplet. On met a jour egalement les parametres de la droite optimale. Avant la mise en oeuvre de 
I'etape A, il faut initialiser la valeur Eop a une valeur faibte negative, par exemple -1. 

Dans ce cas. la rapidite d'obtention de la droite de regression ne depend pas du mode de scrutation 
des triplets. Sa complexite est de degre Une mise en oeuvre materielle peut beneficier de la grande 
regularite de I'algorithme mis en oeuvre. It est peu sensible aux troncatures des donnees et fournit une 
45 solution exacte. 



Troisieme variante de la premiere phase de la methode. 

Dans cette troisieme variante (figure 5), on selectionne d'abord une paire de points a laquelle on ajoute 
50 un point supplementaire, situe entre ces deux points, afin de former un triplet de points. Pour cela. on 
modifie les etapes A. B et C de ta premiere vahante, les autres etapes restant les memes. Les etapes 
modifiees sont teltes que : 

Al - modifie I'etape A en operant une selection de deux couples de valeurs appartenant a ladite 
suite, tel qu'il existe au moins un couple additionnel intermediaire ayant une grandeur indepen- 
55 dante <Xi - Xs) comprise entre les grandeurs independantes dudit couple pour constituer au 

moins un triplet de couples. 

A2 - modifie I'etape A premierement en determinant une fonction lineaire annexe qui contient les 
deux couples selectionnes et deuxiemement en determinant des secondes erreurs entre les 
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granaeurs dependantes des couples intermediaires possibles et ladite fonction lineaire annexe : 
et lorsque ces secondes erreurs sont toutes de meme signe. selection du couple 
internnediaire fournissant la plus grande seconde erreur, en valeur absolue. pour fornner un 
triplet de couples de valeurs forme du couple internnediaire et des deux couples 
5 selectionnes. 

. et lorsque ces secondes erreurs sont de signes differents, reprise de la nnethode a Tetape 
A1, 

Bi • I'etape B est effectuee avec ledit triplet selectionne. 

Cl - modifie I'etape C en selectionnant un couple additionnel dont la grandeur independante n'est 
w pas comprise entre la grandeur independante des deux couples selectionnes. 

Lorsque I'erreur | Eprn | est superieure a l*erreur Et (bloc 110), la methode reprend a I'etape Al (lien 
101) avec une nouvelle selection de paire de points (bloc 102a). 

On observe que le deroulement de cette troisieme variante depend de la scrutation des vateurs et done 
des valeurs elles-memes. La complexite de la mise en oeuvre de cette variante varie entre N et N- ce qui 
75 lui confere un certain avantage par rapport aux variantes precedentes. Une forte parallelisation des moyens 
de mise en oeuvre peut etre operee mais I'implementation des moyens peut presenter un manque de 
regularite ce qui peut constituer un handicap pour une realisation integree. Cette variante est peu sensible 
aux erreurs d'arrondis des valeurs et fournit une solution exacte. 

20 Quatrieme variante de la premiere phase de la methode. 

Elle concerne la determination de la droite de regression a partir des enveloppes. 

II est possible de reduire le nombre de triplets a examiner en determinant des enveloppes respective- 
ment ^uperieure et inferieure entourant les points extremes dans la representation bidimensionnelle de 
25 I'ensemble de points. Une enveloppe superieure ou une enveloppe inferieure est definie telle qu*en joignant 
par une droite deux points adjacents quelconques de Tenveloppe, tous les autres points soient situes d'un 
meme cote respectivement de I'enveloppe superieure ou de I'enveloppe inferieure. On determine ainsi tous 
les points appartenant a ces dites enveloppes. 

La determination de la droite de regression va consister a considerer les paires de points adjacents 
30 d'une des enveloppes auxquels on associe un point intermediaire n'appartenant pas a ladite enveloppe 
pour constituer un triplet et operer comme cela vient d'etre decrit dans le cas des paires de points de la 
troisieme variante. Si une solution optimale n'a pas ete trouvee. on considere les paires de points de Tautre 
enveloppe. 

Pour mettre en oeuvre une enveloppe. on selectionne une paire de points adjacents appartenant a 
35 I'enveloppe. On determine alors s'ii existe un point dispose de telle fagon que son abscisse soil 
intermediaire entre les abscisses des points selectionnes. Lorsque ce point n'existe pas on passe a une 
autre paire de points de ta meme enveloppe. Pour certaines paires, lorsqu'il apparail qu'il existe un ou 
plusleurs de ces points Intermediaires. on cholsit le point intermediaire le plus eloigne de la droite 
contenant les deux points de ta paire pour former un triplet et pour determiner une droite de regression. 
40 Pour determiner si cette droite de regression peut etre selectlonnee comme droite de regression optimale, 
ta methode met en oeuvre les memes operations que celles decrites prealablement dans le cas de la 
troisieme variante. 

Pour cela on modifie la troisieme variante telle que (figure 6), prealablement a I'etape 102a (figure 5), la 
premiere phase de la methode comprend une etape (bloc 100a) de determination d'une enveloppe 

45 inferieure et ou d'une enveloppe superieure reunissant les points les plus extremes de I'ensemble de 
points, la selection des paires de points a I'etape 102a etant faite parmi les points adjacents appartenant a 
Tune ou i'autre enveloppe. La selection de ladite paire de points est effectuee lorsqu'il existe au molns un 
point intermediaire ayant une abscisse situee entre les abscisses des points de la paire de points. S'il 
existe plusleurs points Intermediaires. on forme le triplet avec le point intermediaire le plus eloigne de la 

50 droite passant par les deux points qui forment ta paire de points. Si une solution n'est pas trouvee avec la 
premiere enveloppe, on poursuit le traitement avec la seconde enveloppe. 

La determination des enveloppes Inferieure et superieure est effectuee selon I'organlgramme de la 
figure 7. L'Indice des points allant croissant avec la variable x, le premier point Po fait ainsi partle des deux 
enveloppes. Considerons d'abord I'enveloppe inferieure, les points appartenant a I'enveloppe Inferieure sont 

55 reperes par la lettre Q. Un point courant Qv est repere par I'indlce v. 

A I'etape 400 les deux premiers points de I'enveloppe sont : Go = Po et Qi = P, . Le dernier point 
courant Qv est repere : v = i pour Pi . Le nombre k de points traites est comptabillse. 
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Etape 402 : un premier test est effectue pour determiner si le dernier point Pn-i a ete traite pour 
detecter la fin de la determination de Tenveloppe. 

Etape 404 : dans te cas contraire on teste si v > i . Lorsque v < l on incremente v (v = v + l ) et on 
prend le point Pr comme point Qv (etape 407). On incremente I'indice k pour traiter le point P suivant 
5 (etape 409) at reprise de la methode a I'etape 402. 

Lorsque v ^ 1. on calcuie una droite passant par les points Qv - i at (etape 406) et on determine le 
signe « de Terreur sur la variable dependante entre le point P^ courant et cette droite (Qv -i, Qv) Ceci a 
pour but de determiner si te point courant est au-dessus ou en dessous de la droite (Qv - i. Qv). 

Lorsqua le signe € S o. il faut supprimer le dernier point Qv. decrementer v tel que v = v - 1 (etape 
w 410) et reprendre la methode a I'etape 404. On peut de cette maniere etre conduit a supprimar certains 
points deja acceptes lorsqu'un point suivant necessite de les supprimer. 

Lorsque le signe € est strictament positif. la methode reprend a I'etape 407 avec un point suivant. 
Get organigramme s'applique a la determination des enveloppes inferieure et superieure en inversant le 
signe de I'erreur a considerar. 
15 Pour faire comprendre les mecanismes ainsi mis en oeuvre considerons, a titre d'exemple. le cas 
simple forme par les points Pi. P2. P3. P4 de la figure 2 et detarminons I'enveloppe inferieure. Le point P: 
est le premier point de I'enveloppe d'ou Qo = Pi . On prend P: comme deuxieme point Qi . On calcuie la 
droite Lo reliant Pi, P2. Le segment Pi, P? est accepte comme premier segment de I'enveloppe basse. On 
examine si P3 est au-dessus de la droite Lo. ce qui est le cas sur la figure 2. Le segment P2 P3 est accepte 
20 comme deuxieme segment de I'enveloppe basse. On calcuie ensutte la droite Li reliant Pj et P3 et on 
examine si P* est au-dessus de la droite Li . Ce n'est pas le cas sur la figure 2. Ceci entraine que le point 
P3 et la droite Li sont d'abord ecartes. De meme, P* est en dessous de la droite Lo done le point P2 et la 
droite Lo sont ensuite aussi ecartes comme ne pouvant pas appartenir a I'enveloppe inferieure et les 
operations se poursuivent avec la determination de la droite L? reliant Pi. P*. puis I'examen des points 
25 sutvants (non representes). 

La complexite de mise en oeuvre de cette variante basee sur des enveloppes est de degre et 
depend de t'ordre de scrutation des donnees. Cette complexite est moindre que celle des variantas 
precedentes et de ce fait fournit un resultat rapidement. La regularite de I'implementation est moyenne mats 
cette variante est peu sensible aux erraurs d'arrondi des valaurs et delivre une solution exacte. 
30 Pour certaines applications, f peut etra souhaitable d'accroilre la precision de la determination de la 
fonction d'approximation dans certains domaines de la grandeur indepandante x et d'affecter des coeffi- 
cients da ponderation Wj aux points P;. C'est par exemple le cas lorsque la fonction d'approximation est 
faiblement variable avec la grandeur indepandante x. A certains points peuvent alors etre affectes das 
coefficients de ponderation. Ceux-ci peuvent etre communs a plusieurs points ou etre individuets pour 
35 chaque point. Cas coefficients de ponderation Wi sont, par la suite, consideres strictement positifs. 

Dans les cas ou des coefficients de ponderation existent, on definit une erreur Epo entre un point P, et 
la droite de regression D tel que : 

EpoiP.D) = lV,[y, . (p.X, + q)] 

40 

ou EpD est une valeur signee. La determination de la droite de regression D pour trots points P„ P^, P^ est 
alors modifiee en ce que les quantites p et q de I'equation (I) daviennent : 

(2) p = NUMP'DET et q = NUMQ/DET 

ou les quantites NUMP. NUfVlQ et DET sont definies par : 

DET = Wi.W,.(Xi.X,) + Wi.WK.(X,-Xj) + W..W,.(Xk-X,) 
NUMP = Wi.Wi.(Y,-Y,) + Wj.W,.(Y,-Yj) + Wi.WK.(Y,-Yi) 
50 NUMQ = W..Wj.(Xj.YrX;.Yj) + W,W,.(X,.Yj-Xi.YO + W..W,.(X,.Y.-X..Y,) 

Par ailleurs. I'erreur Ej associee a ce triplet peut etra exprimee et calculee par : 
avec 

55 Ej{P.P„P,) = \Epo{P,.D\ 

EpD{Pi,D) = -EpdP„D) = Epd{P,.D) 

ou 
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DET 

5 

La premiere et la seconde variante de la methode de la premiere phase decrites precedemment 
peuvent etre mises en oeuvre en faisant intervenir les coefficients de ponderation ci-dessus. Cette mise en 
oeuvr^ peut etre effectuee en programmant un calculateur. 

Dans le cas ou Ton selectionne une paire de points en faisant intervenir des coefficients de ponderation 
70 W, la methode presente les amenagements suivants. 

Pour les points Pj, P|, Py qui forment le triplet, on definit un facteur F, qui lie les coordonnees 
respectives et les coefficients de ponderation respectifs des points P,, P|, Pk. Au point central 1, on associe 
le facteur F, tel que : 




20 Ce facteur F, intervient dans la determination du point intermediaire P| qui sera a selectionner (i < I < k) 
pour constituer le triplet en tenant compte des poids affectant les points. On choisit un point Pi et on 
calcule la droite de regression Di associee au triplet P*. Pi. P^ et Terreur Eti associee au triplet. 

Cette variante de la methode consiste a modifier uniquement I'etape A2 de la troisieme variante (bloc 
102b, figure 5). Cette etape determine I'existence et la valeur d'un point intermediaire servant a former un 
25 triplet. On cherche d'abord a former une droite de regression situee en dessous de Pj et de P^. Pour 
chaque point intermediaire P|, on determine si F, = 1 et Epo (Pi. D,) < 0. Si au moins un point verifie cette 
condition, il n'existe pas de droite de regression situee en dessous de Pj et P^, sinon on determine une 
grandeur Gmax pui est la quantite maximale parmi ; 

. d*une part les quantites Epo (Pi, Dp) pour tous les points intermediaires, 
30 d'autre part les quantites (a) suivantes, uniquement pour les points intermediaires pour lesquels F| < 1 

avec 



35 YT^'^^^ ^^I'^i^ ' 

S'il existe au moins un point intermediaire tel que F, > 1. on definit egalement une grandeur G^jn qui 
est la valeur minimale prise par (a) uniquement pour les points intermediaires tels que F, > 1. On examine 
40 alors s'il existe au moins un point intermediaire tel que : 
Epo (Pi, D,) 2: 0 
Gmax ^ Of Eti 

et Gmin- ot ^ Eti (a, coefficient Z 1). S'il existe un tel point, il est choisi comme point intermediaire pour 
former le triplet. 

45 S'il n'existe pas de point intermediaire tel que Fi > 1, on teste s'il existe au moins un point 
intermediaire tel que : 
EpD (Pi. D.) ^ 0 

Gmax ^ ci Eti (o : Coefficient i l). 

Si un tel point existe, il est choisi comme point intermediaire pour former le triplet. 
50 Si aucun triplet n'a ete constitue. on cherche a former une droite de regression situee au-dessus des 
points Pj et P^. La meme methode est reprise en inversant le signe des erreurs Epo- 

Si aucun point P| ne peut etre selectionne, on recommence avec une autre paire Pj. Pn. 
Lorsque I'ensemble de points a traiter est trop important pour etre represente par une seule droite de 
regression, la methode code alors plusieurs droites de regression chacune etant determinee selon la 
55 methode decrite precedemment. 

La figure 13 represente un exemple dans lequel ia fonction d'approximation est constituee de plusieurs 
droites de regression. 
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Dans une premiere situation, de par la connaissance des grandeurs a traiter. on peut vouloir loiposer 
des linnites a chaque droite.de regression suivant la grandeur independante x. Amst on peut vouloir 
disposer d'une droite de regression Da entre les valeurs [X^.X^i de la grandeur, borne Xa incluse. borne Xt, 
exclue. De meme avec Dt, pour [X^,, Xc[ et D^ pour [Xc. X^^ . Dans ce cas le probleme revient a deiernniner 
5 une droite dans un donnaine limite et a appliquer a chaque fois la methode deja decnte. 

Mais il est possible de faire que pour chaque droite, la methode determine des iimites suivant ta 
grandeur independante x sans qu'elles soient fixees au depart. 

Le principe de la determination d'une limite optimale entre deux drouv^s de regression adjacentes est 
represente sur la figure 14-A. Soient deux droites de regression Dl et D2 non optimales. La droite D1 est 
70 determinee a partir de Nl points et la droite D2 est determinee a partir des N2 points restants avec N1 + 
N2 = N. ou N est le nombre total de points. Les droites Dl et D2 ferment une approximation avec 
respectivement des erreurs maximales El et E2. Supposons que E2 soit superieure a El et que Ton 
cherche a reduire E2, il faut alors transferer vers Dl un point appartenant a D2. Lorsque le nombre de 
points diminue. Terreur resultante est soit inchangee soit reduite. Lorsque te nombre de points augmenie. 
15 I'erreur resultante est soit inchangee soit augmentee. En fonction du nombre n de points appartenant aux 
droites D1 et D2. les variations des erreurs El et E2 peuvent presenter les allures representees sur la figure 
14-B pour Texemple choisi. II s'ensuit que I'approximation obtenue gfobatement par D1 et D2 est optimale 
lorsque ces erreurs sont voisines, c'est-a-dire situees dans la zone delimitee par un trait discontinu sur la 
figure 14-B. Ces erreurs ne sont pas necessairement egales compte tenu du caractere discontinu des 
20 valeurs de la grandeur independante. 

Pour determiner la valeur d'abscisse Xii^ entre les deux droites : 

- on determine la droite Dl sur un certain nombre de points et on calcule I'erreur El maximale, 

- on determine la droite D2 sur les points restants et on calcule I'erreur E2 maximale, 

- on compare El et E2 et on transfere un point de la droite qui presente la plus forte erreur vers ta 
25 droite qui presente la plus faible erreur. 

- on determine la valeur limite X^^ lorsqu'il se produit une inversion dans le rapport entre les deux 
erreurs. 

Cette methode s'applique par recurrence a plusieurs droites delimltees entre elles par plusieurs Iimites 
d'abscisses. 

30 Les droites sont determinees a partir d'une suite discrete et limitee de grandeurs de mesure. 
Neanmoins. pour I'exploitation des droites de regression, il est necessaire de definir leur domaine 
d'existence qui s'etend sur un continum de valeurs situees entre deux Iimites d'abscisses. Or la determina- 
tion des droites fournit une suite de droites qui ne se raccordent pas necessairement par leurs extremites. II 
peut apparaitre utile pour certaines applications d'eviter qu'un ecart apparaisse sur la valeur de la grandeur 

35 dependante y pour des valeurs voisines (X„m - f) et {Xum + «) de la grandeur independante ou t est une 
tres faible valeur. II est possible de choisir que I'abscisse limite X),„, appartienne exclusivement a Tune ou a 
I'autre droite. II est aussi possible de faire suivre la determination des suites de droite d'une procedure de 
raccordement de droites. 

Ceci est represente sur la figure 15. Une solution peut consister a valider la droite Dl jusqu'a I'abscisse 
40 du premier point appartenant a la droite D2 et a recalculer la droite D2 a partir de la nouvelle valeur de 
grandeur dependante y ainsi determinee. En conservant le dernier point appartenant a D2. on determine 
ainsi une nouvelle droite D'2 representee en pointille sur la figure 15. Une procedure analogue peut etre 
appliquee pour substituer la droite D'3 a la droite D3, On obtient ainsi un ensemble forme de plusieurs 
droites de regression forme, dans I'exemple decrit, par les droites Dl. D'2, D'3. Get ensemble constitue une 
45 approximation des grandeurs de mesure en reduisant une erreur maximale entre les grandeurs de mesure 
et Tensemble de droites. Une variante plus adaptee a cette reduction consiste a faire que ce soient les 
droites correspondant aux plus fortes erreurs qui imposent leurs points Iimites comme points Iimites aux 
droites correspondant a des erreurs plus faibles. 

50 Dispositif cour une mise en oeuvre de la premiere phase 

Un schema d'un dispositif pour generer une fonction d'approximation selon I'invention est represente 
sur la figure 8. Le dispositif 5 regoit des couples (X^, Y,) associant la grandeur dependante Y, a la grandeur 
independante X,. Les couples (X;. Y,) entrent dans les premiers moyens 10 pour determiner et coder la 
55 fonction lineaire de regression formant une approximation des couples de mesure. Les codes specifiques 
ainsi determines sont transmis (connexion 9) aux seconds moyens 17 qui determinent des seconds couples 
(Xa. Y'a) d'apres les codes specifiques a la fonction d'approximation. 
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On distingue d'une part la seconde variante qui determine toutes tes droites de regression passant par 
chaque connbinaison realisable de triplets avec la suite de premiers couples, et les autres varianies pour 
iesqueltes a chaque droite de regression determmee (associee a son erreur Et), on examine si les points 
additionnels restanis fournissent bien des erreurs additionnelles inferieures ou egales a Terreur de triplet Et- 
5 La figure 9 represente un dispositif adapte a la seconde variante. II comprend : 

- une memoire 12c MEM qui stocke notamment tous les points appartenant a I'ensemble de points a 
traiter. Ces points sont representes par leurs coordonnees (x, y). et eventuellement leurs poids W ou 
leurs inverses lAA/, 

• une unite de calcul 13c COMPUT qui calcule, pour chaque triplet selectionne. la droite de regression 
10 adaptee a chaque triplet, c'est-a-dire les codes p. q de la droite et Terreur Et associee au triplet. 

De plus, un controleur lie CONTR permet d'effectuer ta gestion des operations et d'adresser des 
nouveaux triplets en effectuant les lectures-ecritures de la memoire 12c et le chargement de I'unite de 
calcul 13c par de nouveaux triplets. La selection de la droite de regression qui est a conserver, c'est-a-dire 
pour cette variante celle delivrant la plus grande erreur de triplet, est effectuee par Tunite de calcul I3c. 
15 La figure 10 correspond au cas des autres variantes pour lesquelles pour chaque droite courante de 
regression on examine si les autres points additionnels de la suite delivrent une erreur inferieure a I'erreur 
de triplet. 

Pour cela. les premiers moyens 10 comprennent une unite de calcul 13c et une unite de comparaison 
14c formant des moyens de calcul I9c. 
20 L'unite de calcul 13c, transmet les codes p. q de la droite de regression du triplet courant a I'unite de 
comparaison 14c COMPAR qui determine si les points additionnels de I'ensemble de points generent une 
erreur plus faible avec cette droite de regression que celle generee par les points du triplet courant. Pour 
cela. I'unite de comparaison I4c effectue le test : 

25 I Ep^ |'> Et 

Si la droite de regression courante n'est pas acceptee {test positif). un autre triplet est selectionne et une 
procedure analogue est a nouveau effectuee. Si le test est negatif pour tous les points additionnels. la 
droite de regression est acceptee et ses parametres sont charges par I'unite de calcu! 13c dans la memoire 
30 12c. 

Pour les variantes qui operent a partir d'une paire de points. I'unite de calcul 13c est programmee pour 
determiner les points intermediaires et pour former les triplets issus des paires de points, ceux-ci etant eux- 
memes eventuellement issus d'enveloppes de points. Dans ce dernier cas. I'unite de calcul 13c est 
egalement programmee pour determiner les enveloppes. L'unite de comparaison I4c effectue ensuite les 
35 comparaisons des erreurs se rapportant aux differents triplets. 

Pour effectuer ta gestion des operations et adresser des nouveaux triplets (lorsque des triplets 
examines n'ont pas donne satisfaction), le dispositif comprend. comme precedemment. un controleur lie 
qui gere : 

- les lectures/ecritures de la memoire I2c, 

40 - le chargement de I'unite de calcul 13c par de nouveaux triplets. 

- I'appel vers I'unite de comparaison I4c de tous les points additionnels a examiner successivement. 

II est possible de mettre en oeuvre la premiere phase de la methode a I'aide de moyens de 
comparaison 14c (figure 12) ayant une structure neuronale. Considerons le cas te plus complet ou il existe 
des coefficients de ponderation associes a chaque point Pj. Lorsque ces coefficients n'existent pas il suffit 
45 de teur donner une valeur unite dans les explications qui vont suivre. A litre d'exemple. on peut stocker 
dans la memoire 12c (figures 9. 10) les inverses des coefficients Wj affectes individuellement aux points Pj. 
L'unite de calcul 13c determine selon les equations (2) et (3), pour un triplet donne, les codes : -p. -q Et et 
-Et- 

La condition a tester dans I'unite 14c pour un point additionnel P^ a tester, de parametres X„„ Y^. W^. 

50 est : 

Y„, - pX^ - (Et-W^) - q > 0 
ou 

55 

- p.X^ + (Et/W^) - q < 0. 
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Ces tests sont aisement reahsables par une unite I4c ayant une organisation neuronals. En effet. ce 
test est constitue de fonctions lineaires. de fonctions de seuillage et de fonctions logiques aisement 
realisees dans un reseau de neurones. 

Une telle unite neuronale de comparaison est representee sur la figure 12. Elle comprend trois 

5 neurones N1, N2. N3. Les neurones Ni et N2 regoivent des donnees El. E2. E3. E4. Le neurone N3 regoit 
les sorties des neurones NI et N2. Chacune des entrees de ces neurones est affectee d*un coefficient 
synaptique Ci selon la technique connue mise en oeuvre dans les reseaux de neurones. Cette technique 
est par exemple decrite dans R.P. LIFPMANN. "An introduction to computing with neural nets" IEEE ASSP 
Magazine, April 1987, pp. 4 to 22. 

10 Pour effectuer les tests precedents, on programme les neurones NI et N2 selon le tableau I : 

TABLEAU I 



15 



Neurone 


E, 




E2 




E3 


C3 


E4 


C* 


1 




1 




- P 




-Et 


1 


-q 


2 




1 




- P 


1>W^ 


Et 


1 


- q 



20 Tous les coefficients synaptiques du neurone N3 sont egaux a 1 . 

Les neurones Nl et N2 calculent chacun un potentiel neuronal s tel que : s = C,.E, 
A ce potentiel s, il faut appliquer une fonction deactivation A. Pour les neurones NI et N3, la fonction 
d'activation Ai est : 

25 At (s) = 1 si s > 0 
Ai(s) = 0 si s S 0. 

Pour le neurone N2, la fonction d'activation A2 est : 

30 A2(s) = 1 si s < 0 

A2(s) = 0 si s ^ 0. 

La sortie du neurone N3 vaut 1 si la condition a tester est satisfaite et vaut 0 dans le cas inverse. 

35 On observe sur le tableau I que les donnees -p, -q. -Et, + Et sont des codes qui apparaissent comme 
coefficients synaptiques dans le schema de la figure 12. Les donnees Y^. X^, I'W^ sont des donnees qui 
apparaissent sur les entrees du meme schema. Lorsque tous les points addittonnels ont ete testes et 
qu'une droite de regression a ete selectionnee les codes p. q. W charges dans la memoire I2c sont alors 
utilises pour mettre en oeuvre la seconde phase de la methode (figure 11). 

40 L'interet presents par la realisation neuronale decrite est de pouvoir paralletiser les differentes 
operations a mettre en oeuvre selon les vahantes deja decrites. Le fonctionnement d'une telle realisation 
neuronale est alors tres rapide. 

Lorsque les droites de regression ont ete determinees. la seconde phase consistant a calculer des 
seconds couples (Xa, Y'a) de valeurs des grandeurs est alors mise en oeuvre dans des moyens de 

45 decodage 17 (figure 11). Les codes des droites sont charges dans une memoire 12a qui au cours de la 
seconde phase est adressee. par un controleur 11a, pour fournir les codes des droites de regression 
adressees. La memoire 12a, organisee en lignes par exemple. contient pour chaque droite de regression 
les parametres p, q et xl oCj xl est la limite suF>erieure d'abscisse pour laquelle chaque droite de regression 
est definie. 

50 La memoire 12a contient ainsi un tableau de parametres 



^Lm* Pm-l- 

correspondent aux m droites de regression stockees. 

12 
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Les moyens de decodage 1 7 comprennent : 

- le controleur n^. 

- la memoire I2a, 

- et une unite de decodage I3a. Celle-ci regoit une valeur de requete Xa de la grandeur independante 
5 pour laquelte on desire obtenir un resultat Y'a issu de la fonction d'approxinnation. Pour cela le 

controleur 1 1^ presente successivement certaines lignes de codes de la nriemoire I2a et, pour chaque 
ligne de codes. Tunite de decodage i3o teste le code xi pour determiner si la valeur Xa de la 
grandeur d'entree peut etre traitee par )a droite de regression codee par ladite ligne. Si ceta n'est pas 
le cas, I'unite 13a requiert du controleur 11 a une autre ligne de codes. Lorsque cette condition est 
w verifiee, I'unite 13a e^fectue le calcul : 

Y'a = Po-Xa + qn 

si la ligne d'ordre n est la ligne adequate. 
15 Connme cela vient d'etre decrit, les premiers moyens 10 peuvent constituer un dispositif de codage et 
les seconds moyens 17 peuvent constituer un dispositif de decodage, les deux dispositi*s etant separes. lis 
peuvent en effet etre distants I'un de I'autre, les operations de codage etant separees des operations de 
decodage. Mais les deux moyens 10 et 17 peuvent aussi former un dispositif de codage/decodage unique. 
Dans ce cas certains organes peuvent etre affectes successivement a plusieurs operations. Notamment : 
20 - il peut s*agir d'une memoire 12 unique qui peut stocker les parametres des couples de valeurs et les 
codes des droites de regression determinees. Les unites 12a et 12c torment alors la memoire 12, 

- II peut s'agir d'un controleur 1 1 unique qui regroupe les fonctions des controleurs 1 lc et 1 la. 

- I'unite de calcul I3c et I'unite de decodage I3a peuvent egalement former une seule unite 13. 

II est egalement possible de realiser les moyens de decodage 17 (figure 11) selon une organisation a 

25 structure parallele. Une solution consiste alors a comparer, en parallele. la grandeur de requete Xa h tous 
les codes xl. On obtient un ensemble de signaux Tl representant les resultats de ces tests de comparaison. 
On realise ensuite. en parallele, pour toutes les valeurs de L, un test qui combine deux a deux les signaux 
Tl et et qui indique si la grandeur de requete Xa est dans I'intervalle d'indice L. Ce test est verifie pour 
une unique valeur de L. notee n. On accede alors aux codes pn et qn correspondants. et on en deduit la 

30 grandeur Y'a de la meme maniere que precedemment. 

Les premiers moyens 10 (de codage) associes aux seconds moyens 17 (de decodage) peuvent etre 
utilises pour determiner la valeur d'une fonction d'approximation par au moins une droite de regression. 
Ces determinations peuvent etre effectuees pour des valeurs quelconques de la grandeur independante (a 
I'interieur de bornes predefinies definissant I'etendue de Taction de chaque droite de regression). Cette 

35 methode evite de stocker inutilement des tables de valeurs, pour lesquelles toutes les valeurs ne seront pas 
exploitees. Selon la methode on ne determine que les seules valeurs necessaires a I'application. La 
methode selon ['invention presente I'interet de ne calculer que les valeurs necessaires. Le dispositif de 
invention peut etre utilise en combinaison avec un processeur neuronal pour calculer une approximation 
d'une fonction non lineaire par exemple une fonction sigmoVde. Un processeur neuronal calcule des 

40 potentiels neuronaux issus de sommes de produits de coefficients synaptiques par des etats de neurones. 
La technique est developpee dans le document de R.P. LIPPMANN deja cite. Chaque potentiel neuronal 
doit etre soumis a Taction d'une fonction non lineaire. Pour cela. chaque potentiel neuronal constitue la 
grandeur independante decrite precedemment et le calcul de la fonction th(x) par exemple est effectue pour 
les seules valeurs de requetes utiles. La methode fournit des valeurs approximees. mais cette approxima- 

45 tion peut etre obtenue avec une precision qui peut etre accrue selon le nombre de droites de regression 
calculees et eventuellement les coefficients de ponderation W. La precision reste neanmoins tributaire de la 
precision des coordonnees des couples de depart. 

Une telle methode est particulierement interessante pour les calculs de fonctions connues (telles que 
des fonctions mathematiques) ou de fonctions inconnues explicitement. par exemple une fonction represen- 

50 tee par des points de mesures. que Ton desire simplifier par une fonction lineaire de regression. 

Dans sa version neuronale. I'invention est interessante dans des applications neuronales car elle 
apporte non seulement un traitement homogene mais aussi une grande compacite de I'architecture 
materielle necessaire. 

II est possible de combiner les variantes de la methode. dans leurs mises en oeuvre materielles pour 
55 combiner leurs avantages. Ainsi par exemple, les combinaisons suivantes sont avantageuses : 
Dans le cas ou les poids existent, deux combinaisons apparaissent interessantes : 
Dans la premiere combinaison. on utilise tout d'abord la premiere variante basee sur des triplets de 
points. Si une solution n'est pas trouvee. on enchaTne avec la deuxieme variante basee aussi sur des 
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triplets de points. Ceci permet d'obtenir les avantages suivants : 
Queues que sotent les donnees. on oblient une solution. 

Dans le cas oO on exige qu'on obtienne toujours une solution exacte, cette connbinaison pernnet 
d'obtenir la rapidite nnaximale. En effet. dans le cas general la solution est fournie dtrectennent par 
5 la prenniere variante appelee, qui est cede fournissant le plus rapidement une solution exacte dans 

le cas ou les poids existent. La deuxieme variante est beaucoup plus lente mais elle permet de 
garantir que la combinaison fournit toujours une solution. 
Dans la deuxieme combinais::n. on utilise tout d'abord la variante basee sur des paires de points. Si 
elle ne trouve pas une solution (qui peut etre approchee). on utilise la deuxieme variante basee sur 
10 des triplets de points. Ceci permet d'obtenir les avantages suivants : 

Quelles que soient les donnees, on obtient une solution (eventuellement approchee). 
Cette combinaison est en moyenne plus rapide que la precedente. 
En contrepartie. elle ne fournit generalement qu'une solution approchee. 
Dans le cas ou les poids n'existent pas, une combinaison interessante consiste a utiliser tout d'abord la 
75 variante basee sur Tenveloppe des points. Si elle ne trouve pas de solution, on utilise la deuxieme variante 
basee sur des triplets de points. Ceci permet d'obtenir les avantages suivants : 
Quelles que soient les donnees. on obtient une solution, 
Cette combinaison permet d'obtenir la meilleure rapidite moyenne. 
Lorsque la fonction lineaire de regression a ete determinee selon la methode decrite. il est possible de 
20 determiner non seufement des grandeurs y\ ressortissant a ladite fonction mais egalement des grandeurs 
ressortissant a la derivee de ladite fonction. En effet, pour une valeur de requete Xa de la grandeur 
independante, il suffit d'extraire des moyens de memorisation la valeur p ressortissant a la valeur Xa pour 
connaTlre la valeur de cette derivee correspondante. 

25 Revendications 



1. Dispositif (5) pour generer une fonction d'approximation fondee sur des premiers couples ((Xi. Yi) a 
(Xe, Yg)) de valeurs associant une grandeur dependante (Yi a Y^) a une grandeur independante (Xi a 
Xg), et pour determiner des seconds couples (Xa, Y'a) de valeurs desdites grandeurs d'apres ladite 
30 fonction d'approximation caracterise en ce que le dispositif comprend : 

- des premiers moyens (10) : 

- pour determiner iterativemant au moins une fonction lineaire courante de regression en 
rendant egales. en valeur absolue. des premieres erreurs, de signes alternes. mesurees entre 
des premieres valeurs (Y3, Yc, Y5) de la grandeur dependante pour trois couples (X3. Y3) (X*. 

35 Y4) (X5, Ys) d'une suite desdits premiers couples, et respectivement des secondes valeurs 

(Y*3, Y'a, Y'=) de la grandeur dependante determinees, d'apres ladite fonction lineaire courante, 
pour les memes valeurs (X3. Xa, X5) de la grandeur independante. 

- pour setectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales, 

40 - et pour coder, a I'aide de codes specifiques (p. q), ta fonction lineaire de regression 

selectionnee, 

- et des seconds moyens (17) pour determiner lesdits seconds couples (Xa. Y'a) a I'aide desdits 
codes specifiques. 

45 2. Dispositif selon la revendication 1 caracterise en ce que les premiers moyens (10) comprennent : 

- des moyens (I2c) pour memoriser les grandeurs des premiers couples et lesdits codes specifi- 
ques, 

- des moyens (I3c) (I9c) pour calculer successivement les codes specifiques (p, q) de fonctions 
lineaires courantes de regression et leurs premieres erreurs, egalisees en valeur absolue, 

50 respectives (Ej) et pour selectionner au moins une fonction lineaire courante de regression. 

- et des moyens (lie) pour controler les moyens (I2c) pour memoriser et tes moyens (I3c) pour 
calculer. 



3. Dispositif (5) selon les revendications 1 ou 2 caracterise en ce que pour setectionner ladite fonction 
55 lineaire, les premiers moyens comprennent des moyens (14c) : 

- pour determiner pareillement des erreurs additionnelles pour les autres couples (Xi, Yi) (X2. Y2) 
(Xs . Ye ) de valeurs de ladite suite, 

- pour comparer les erreurs additionnelles avec lesdites premieres erreurs egalisees (Et), 
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- el pour selectionner la fonction lineaire courante qui fournit des erreurs additionnelles mferieures 
ou egales en valeur absolue auxdites premieres erreurs egalisees (Et). 

Dispositif selon ta revendication 2 caracterise en ce que. pour selectionner la fonction lineaire de 
regression, les moyens pour calculer (I3c) selectionnent cells des fonctions lineaires courantes de 
regression qui delivre des premieres erreurs egalisees (Et) maximales. 

Dispositif selon une des revendications 1 a 4 caracterise en ce que les seconds moyens (17) 
comprennent : 

- des moyens (12a) pour memoriser lesdits codes specifiques d'au moins une fonction lineaire de 
regression. 

- des moyens (13a) pour selectionner. sur requete d'une valeur de grandeur independante (Xa). 
une fonction lineaire de regression ressortissant a ladite valeur. et pour calculer. a partir des 
codes specifiques regus des moyens pour memoriser. une valeur de grandeur dependante (Y'a) 
associee a la valeur de grandeur independante (Xa), 

- et des moyens (lla) pour adresser aux moyens pour calculer. et cur leur demande, les codes 
specifiques des fonctions lineaires de regression. 

Dans un dispositif pour generer une fonction d'approximation fondee sur des premiers couples ((Xi , Yi ) 
- (Xs. Ys)) de valeurs associant une grandeur dependante (Yi - Yg) a une grandeur independante (Xi - 
X&). methode pour generer tadite fonction d'approximation a partir des premiers couples (X:, Yi) (Xe . 
Ye) de valeurs, et pour determiner des seconds couples (Xa, Y'a) de valeurs desdites grandeurs 
d'apres ladite fonction d'approximation, caracterisee en ce que, la methode comprend : 
une premiere phase : 

- pour determiner iterativement au moins une fonction lineaire courante de regression en 
rendant egales. en valeur absolue, des premieres erreurs Epo. de signes alternes. mesurees 
entre des premieres valeurs (Y3, Y*, Y5) de ta grandeur dependante pour trois couples (X3. 
Y3), (X4, Y4) (Xs, Ys) d'une suite desdits premiers couples, et respectivement des secondes 
valeurs (Y'3, Y'4, Y's) de la grandeur dependante determinee, d'apres ladite fonction lineaire, 
pour tes memes valeurs (X3, X*, Xs) de la grandeur independante. 

- pour selectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales. 

- et pour coder la fonction lineaire de regression selectionnee a I'aide de codes specifiques. 

- et une seconde phase pour determiner lesdits seconds couples (Xa. Y'a) a I'aide desdits codes 
specifiques. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes 
suivantes : 

A - selection (102) de trois couples (P3, P4. P5) de valeurs parmi ladite suite. 

B - calcul (104) de la fonction lineaire courante de regression D et determination d'une erreur de 

triplet Er=|£'pD| s'y rapportant. 
C - selection (106) d'un couple additionnel (PI). (P2). (P6), 

D - calcul (108) d'une erreur additionnelle Ep^ entre le couple additionnel et ladite fonction, 

E - et lorsque | Ep^ | < E^IIO) pour le couple additionnel. la methode reprend a I'etape C avec 

un couple additionnel suivant. 
F - et lorsque | Ep^ | > Et pour au moins un couple additionnel (110). la methode reprend a 

I'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite de 

couples. 

G - et lorsque | Ep^ | ^ Er pour tous les couples additionnels (112). la fonction lineaire courante 
de regression est codee et stockee (114) en tant que fonction lineaire d'approximation. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes 
suivantes : 

A - selection (102) de trois couples (P3. P4. P5) de valeurs parmi ladite suite, 

B - calcul (104) d'une fonction lineaire courante de regression D et determination d'une erreur de 

triplet Et = \Epd\ s'y rapportant. 
C - comparaison (210) de I'erreur Et avec une erreur optimale Eop initialisee a une valeur 

strictement negative. 
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D - et lorsque Et > foo . n^'se a jour (114) de I'erreur oplimale Eop en remplagant Eop par Et et 

mise a jour des codes d'une fonction lineatre optimale de regression Dop en rennplagant 

ceux-ci par les codes de la fonction lineaire courante (D). 
E - puis retour (1 16) a I'etape A pour seiectionner trois autres couples (102). 
5 F - et torsque tous les triplets de couples de valeurs de ia suite ont ete testes (116). les derniers 

codes de la fonction lineaire optimale Dop constituent les codes de la fonction lineaire 

d'approximation (1 14). 

9. Methode selon la revendication 8 caracterisee en ce que la premiere phase comprend les etapes 
10 modifiees suivantes : 

Al - modifie I'etape A en operant une selection (I02a) de deux couples de valeurs appartenant a 
ladite suite, te! qu'il existe au moins un couple additionnel intermediatre ayant une grandeur 
independante (Xi - Xs) comprise entre les grandeurs 
independantes dudit couple pour constituer au moins un triplet de couples de valeurs. 
;5 A2 - modifie I'etape A premierement en determinant (102b) une fonction lineaire annexe qui 

contient les deux couples selectionnes et deuxiemement en determinant des secondes 
erreurs entre les grandeurs dependantes des couples intermediaires possibles et ladite 
fonction lineaire annexe : 

et lorsque ces secondes erreurs sont toutes de meme signe. selection du couple 
20 intermediaire fournissant la plus grande seconds erreur, en valeur absolue. pour 

former un triplet de couples de valeurs forme du couple intermediaire et des deux 
couples selectionnes. 

et lorsque ces secondes erreurs sont de signes differents, reprise de la metfiode a 
I'etape Ai . 

25 B1 - I'etape B (104) est effectuee avec ledit triplet selectionne. 

Cl - modifie I'etape C (106) en selectionnant (I06a) un couple additionnel dont la grandeur 
independante n'est pas comprise entre la grandeur independante des deux couples selec- 
tionnes- 

30 10. (Vletfiode selon la revendication 9 modifiee en ce que. prealablement a I'etape Al. lesdits couples de 
valeurs etant representes par des points (PI - P6) dans un espace a deux dimensions, la premiere 
phase de la methode comprend une etape (100a) de determination d'une enveloppe inferieure et/'ou 
d'une enveloppe superieure reunissant les points les plus extremes de I'ensembie de points, la 
selection des deux couples de valeurs a I'etape Al etant faite parmi les points adjacents appartenant a 

35 Tune ou I'autre enveloppe. 

11. Methode selon une des revendications 6 a 10 caracterisee en ce qu'a chaque valeur, prise dans ta 
suite de valeurs de la grandeur independante. est associee un coefficient specifique de ponderation qui 
pondere lesdites premieres erreurs. 

40 

12. Methode selon une des revendications 6 a 11 caracterisee en ce qu'elle genere une fonction 
d'approximation comprenant une suite de fonctions lineaires de regression raccordees entre elles. 
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